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The pattern of anisotropies in the Cosmic Microwave Background depends upon the masses and 
lifetimes of the three neutrino species. A neutrino species of mass greater than 10 eV with lifetime 
between 10^^ sec and 10^'' sec leaves a very distinct signature (due to the integrated Sachs-Wolfe 
effect): the anisotropies at large angles are predicted to be comparable to those on degree scales. 
Present data exclude such a possibility and hence this region of parameter space. For ~ 30 eV, 
r ~ 10^^ sec, we find an interesting possibility: the Integrated Sachs Wolfe peak produced by the 
decaying neutrino in low-fi models mimics the acoustic peak expected in an f2 = 1 model. 
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Introduction. The possibility that one or more of the 
neutrino species is massive has intrigued cosmologists for 
over thirty years A massive neutrino contributes to 
the energy density in the universe. A mass on the or- 
der of 50 eV is sufficient to push the total density in the 
universe up to the level expected if the universe is ffat 
1^. Masses much larger than this are ruled out cosmo- 
logically, based on measures of the age of the universe. A 
caveat to this argument is that an unstable neutrino can 
decay early enough so that its decay products redshift 
away and no longer contribute too much energy to the 
universe. This caveat relies on the fact that the energy 
density of relativistic particles (such as the decay prod- 
ucts) drops faster than that of non-relativistic particles 
(such as massive neutrinos). 

The caveat that neutrinos with mass much greater 
than 50 eV are allowed cosmologically if they decay 
fast enough is important, for present laboratory lim- 
its for the M (r) neutrino are 170 keV (24 MeV) [|. 
The theoretically-predicted decay modes and lifetime 
of a massive neutrino are very model dependent [Q, 
with shorter lifetimes typically arising in familon models 
wherein the decay products are a Majoron and a lighter 
neutrino. Here we consider primarily such decay modes: 
(j), where is a Majoron and V]^ are any 



of Vr, l^u 



Although our limits will also be applica- 



bleQ to photon-producing decays: vu + there are 

already other restrictive limits on such decays. For the 
decay Vr u^-\- (p, a plausible lower limit can be placed 
on the lifetime from upper limits on the branching ratio 
of r ^ ^ + I 



*Our limits also apply to 3-body decays as long as the decay 
products are relativistic. 



'0.67 MeV\'' 
T > I sec 

m^^ J 



(1) 



In this Letter, we argue that a large part of the 
mass/lifetime parameter space is strongly disfavored by 
recent observations of anisotropies in the Cosmic Mi- 
crowave Background (CMB). Dozens of observations 
have confirmed that (i) anisotropies exist in the CMB 
and (ii) the level of anisotropy is higher on degree scales 
than on the large angular scales probed by the COBE 0| 
satellite. Any viable cosmological model must account 
for these observations. We show here that a cosmological 
model with a decaying neutrino with mass m > 10 eV 
and lifetime in the range lO^'^sec < r < lO^^sec fails to 
produce substantial power on degree scales and therefore 
fails to reproduce the current observations. 

Our argument assumes an otherwise standard Cold 
Dark Matter (sCDM) model, wherein the total density 
total is equal to the critical density (fJ = 1). This density 
consists of the energy associated with one species of neu- 
trino with mass m and lifetime t and its relativistic decay 
products; ordinary baryons with an abundance suggested 
by recent measurements of deuterium (flBh^ — 0.02); 
and the remainder in the form of cold dark matter. For 
most of our discussion, the Hubble constant is fixed at 
Ha = 50kmsec~^ Mpc^^. Finally the CMB anisotropy 
spectrum depends on the primordial spectrum, which we 
assume to be Harrison-Zel'dovich, with spectral index 
n = 1. We also assume the decay products of the mas- 
sive neutrino are sterile, i.e. no photons are produced. 
This model serves as a useful framework within which to 
make our argument. After presenting the argument and 
the subsequent limits on m, t, we discuss the changes in 
these limits if the true parameters differ from those in 
our canonical model. 
Integrated Sachs-Wolfe Effect. 

In most theories of structure formation, anisotropies in 
the CMB today reffect conditions in the universe at the 
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time of last scattering, when electrons and protons com- 
bined to form neutral hydrogen. After that time (when 
the universe was of order 10^'^sec old corresponding to a 
redshift z ~ 1100) photons travelled freely through the 
universe. In particular, on angular scales smaller than a 
degree, the anisotropics reflect the fact that at last scat- 
tering, the combined electron/ baryon/photon fluid was 
undergoing acoustic oscillations In Fourier space, the 
temperature peaked at wavenumber fcp ~ vr/rs = 0.024 
Mpc~^ where is the sound horizon at last scattering. 
If we expand the anisotropy spectrum today in terms of 
Legendre polynomials: 



(2/ + i; 

47r 



GF,(cos0), 
(2) 



then the power 1(1 -|- 1)C; peaks at Ip ~ fcp77o = 200 where 
770 ~ '^HqI is the conformal time today. In sCDM, the 
ratio of the power at Ip (on degree scales) to that at Z ~ 10 
(COBE scales) is of order 5:1, and indeed this ratio is 
consistent with present data. 

While the anisotropy spectrum on small angular scales 
does indeed derive from conditions at last scattering, per- 
turbations on larger angular scales enter the horizon only 
after this epoch and they are therefore sensitive to condi- 
tions at late times. In sCDM, the gravitational potential 
due to the cold dark matter is constant. Coupled with 
the assumption of a Harrison-Zel'dovich spectrum, the 
constant gravitational potentials lead to flat power on 
large scales. Variants of sCDM sometimes predict devia- 
tions from this canonical prediction; in particular, unless 
the universe is dominated by non-relativistic particles af- 
ter last scattering, the gravitational potentials will not in 
general be constant. An example of this is an open uni- 
verse, where curvature begins to dominate at late times. 
If the gravitational potential $ is not constant at late 
times, then the resulting contribution to anisotropy |^ 
on a scale I is 



Ai = -2 drj—ji[k{m - rj)] 



(3) 



where 77* is the epoch of last scattering. The impact of 
a varying gravitational potential as described in Eq. ^ 
is called the Integrated Sachs- Wolfe (ISW) effect. One 
important feature of Eq. ^ is the factor of two in front; 
this is to be compared with the factor of 1/3 in front 
of the ordinary Sachs- Wolfe effect. Even a small decay 
in the potential can have dramatic implications for the 
anisotropics. 

If $ starts to vary at conformal time ry^,, the scales 
most affected are those just entering the horizon at that 
time, with wavenumbers k ~ Tr/yy^. The spherical Bessel 
function in Eq. || tells us that these will be projected onto 
angular scales of order I Tr{r]o/riy — 1). It will be useful 



to rewrite this order of magnitude estimate in terms of 
cosmological time t (instead of conformal time rj) . Since 
7] cx t^^^ when CDM is the dominant component, the 
ISW effect peaks at 



4 X 10 'secx 1/3 



- 1 



(4) 



Equation ^ tells us that any disturbance to a flat, non- 
relativistic matter dominated universe at times lO^'^ < 
t < lO^^sec will show up in an enhancement in the 
anisotropy spectrum leftward of the peak at Ip ^ 200. 
As we now show, models in which a massive neutrino is 
unstable with a lifetime in the above range produce just 
such a disturbance. 

Anisotropies in Massive Neutrino Models. 

The decay products of a massive unstable neutrino are 
typically very light, so decays turn non-relativistic en- 
ergy into relativistic energy. This leads to a decay of the 
gravitational potential: the relativistic particles do not 
clump as easily as their massive parents. The decay in 
the gravitational potential, as we have seen, leads to an 
ISW effect at I given by Eq. ^ with ty now replaced by 
the massive neutrino lifetime r. While Eq. ^ gives the 
location of the ISW enhancement, the amplitude of the 
enhancement is proportional to the ratio of energy in rel- 
ativistic decay products to that in CDM. Figure 1 shows 
an example for a 10 eV decaying neutrino with lifetime 
T = 10^^ sec. There are several features of note here. 
First, the energy in decay products peaks at a ~ 0.02 
which, since a oc t'^^^, does indeed correspond to ty ~ r. 
Second, decay radiation peaks at only 10 — 20% of the 
density of CDM. However, due to the large coefficient in 
front of the ISW effect, even this small contaminant of 
radiation produces a large change in the CMB spectrum. 
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FIG. 1. The densities of the components of the universe as 
a function of scale factor for a neutrino with m = 10 eV and 
T = 10^^ sec. Densities are with respect to the the total den- 
sity and the scale factor is set to one today. Solid line is CDM 
which dominates at late times; long dashed line is radiation 
in the form of photons and two light neutrino species; short 
dashed line is the density of relativistic decay products of the 
massive neutrino; and the dotted line is the massive neutrino. 
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Figure 2 shows the power spectra of five difi^erent mod- 
els, normahzed by the peak at I ~ 200. The canonical 
curve is sCDM, for which the peak {I — 200)/plateau 
{I = 10) ratio is 5.6/1. The other curves are decaying 
neutrino models, each with a mass of 10 eV. As the neu- 
trino hfetime gets longer, the peak due to the ISW effect 
moves out to smaller I, in accord with Eq. ^ Note that 
the ISW peak is quite substantial even though a mass 
this small produces relatively little energy in relativistic 
decay products. 
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FIG. 2. CMB Anisotropy Spectra for Decaying Neutrino 
models. Solid line is standard CDM with three species of 
massless neutrinos. Curves are normalized at the first peak 
at I ~ 200. All curves are for a neutrino with m = 10 eV. 

Larger mass neutrinos produce larger ISW effects. 
Hence m > 10 eV is ruled out if the lifetime is such 
that the ISW effect is to the left of the peak. Figure 
3 shows a contour plot of the peak/plateau ratio in the 
m, T plane. Also shown is the curve corresponding to 
57 = 1, and the bound on r from equation (]^). A rough, 
very conservative cut would be to disallow the region in 
which the ratio is less than one. As seen from Figure 3, 
this cut excludes a large region of parameter space that 
would otherwise be allowed. 

A number of other constraints can be placed on such 
long-lived neutrinos. If such neutrinos or their decay 
products dominate the density, they can yield an age for 
the universe which is below bounds from globular cluster 
ages |l^. In Figure 3, these limits lie near the f2 = 1 
curve. More restrictive constraints from structure for- 
mation are also possible but these are considerably 
more model-dependent than our very simple result. 
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FIG. 3. Contours of the peak {I = 200)/plateau (/ = 10) 
ratio for CMB spectra. sCDM predicts a ratio on the order of 
five; data certainly excludes values less than one. The upper 
right region of parameter space produces SI > 1. The lower 
left hatched region is the bound on ^ + from equation 
(0). 

Model Dependence of the Constraint. The contours 
drawn in Figure 3 depend on the underlying cosmological 
model and associated set of parameters. How does the 
constraint hold up as these vary? Certainly, all of this 
work is predicated on the assumption that the primordial 
perturbations are adiabatic. This assumption at present 
seems well-motivated both theoretically and observation- 
ally Raising the Hubble constant to a value more 
consistent with present data lowers the peak. A higher 
^Ib would help here, but our value of 0.08 is about as high 
as possible, given constraints from Big Bang Nucleosyn- 
thesis We have also neglected reionization, but this 
too lowers the peak/plateau ratio even further. We have 
neglected tensor modes, but these too work in the wrong 
direction, falling off at /> 100. Our model assumes that 
the decay products are sterile. However, our results also 
apply to photon-producing decays which occur after last 
scattering; such decays would also lead to a large ISW ef- 
fect. But in the parameter range of interest, the photons 
produced in such decays would seriously distort the well- 
measured thermal spectrum of the CMB or produce 
an (unobserved) diffuse photon background [T^ . 

One parameter which could serve to produce a rise on 
small scales is the spectral index of the primordial per- 
turbations, n, which we have set to one. Most inflation- 
ary models |l^ predict values slightly smaller than one, 
which would again reduce even further the peak/plateau 
ratio. However, there are some models (e.g. hybrid in- 
flation [Ts] ] and supernatural inflation jl^ ] ) which predict 
"blue" spectra (n > 1). Even if we allow n > 1, though, 
the vast majority of the excluded m, r region in Figure 
3 will remain excluded. To see why, consider Figure 4, 
which shows the spectrum for a set of excluded m, r val- 
ues. The spectrum can be modified to go through the 
data at / 200 if n — 2.15. Such an extreme value of n, 
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however, is strongly disfavored at both higher and lower 
values of I. 




We thank U. Seljak and M. Zaldariagga for the use of 
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NASA Grant NAG 5-7092. 



FIG. 4. Impact of varying the scalar spectral index. The 
solid line is the sCDM prediction. The dashed line is for 
a neutrino with m — 30 eV and r = 10^^ sec, while the 
dashed-dotted line is the same decaying neutrino model, but 
with a scalar spectral index n — 2.15. This value was chosen 
to make the CMB exhibit a rise at Z ~ (100—200), as indicated 
by observations, shown here by the points. Note the poor fit 
of the tilted model at low and high values of I. 

Finally, we have assumed a flat universe today (SI — 1). 
An open universe, however, would only exacerbate this 
effect since the peak/plateau ratio is already relatively 
small in open models. This is true for two reasons. First, 
curvature-induced gravitational potential decays cause 
extra power at I ^ 10. Second, lowering shifts the 
first acoustic peak to the right, lowering the anisotropy 
at / ^ 200. There is one interesting caveat to this ar- 
gument. In an open universe, the ISW effect due to the 
decay products also shifts to larger /. For small enough 
lifetimes, the ISW peak could be shifted all the way to 
I 200, mimicking the first acoustic peak in a flat uni- 
verse! For example, a model with m = 30 eV, t = 10^'^ 
sec and = 0.3 gives an ISW peak near I — 200 and a 
peak/plateau ratio 6, in rough agreement with current 
observations. 

Conclusions. An unstable neutrino with mass greater 
than 10 eV and lifetime lO^'^ ^ ''' ^ 10^^ sec is ruled 
out by current CMB observations. These limits are quite 
robust to changes in the underlying cosmological model. 
Moreover, they apply generally to any massive particle 
species (e.g. a neutralino or gravitino) that contributes 
significantly to the energy density and decays in the post- 
decoupling time frame. They also rule out the possibility 
of a radiation dominated Universe, proposed in 1984 to 
solve the problem |0. Moving beyond current ob- 
servations, precision measures of the CMB spectrum out 
to / ~ 2000, as expected from the MAP and PLANCK 
satellites [|8|, are capable of probing the mass/lifetime 
plane for all 771 > 1 eV ||ll. 
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